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Abstract. We show decay estimates for the propagator of the discrete Schrodinger 
and Klein-Gordon equations in the form \\U(t)f\\i<*> < Ci 1 + l*l) _<i/3 ||/llii- Thi s 
implies a corresponding (restricted) set of Strichartz estimates. Applications of 
the latter include the existence of excitation thresholds for certain regimes of the 
parameters and the decay of small initial data for relevant V norms. The analytical 
decay estimates are corroborated with numerical results. 



A sequence of (time evolving) harmonic oscillators which interact only with their 
immediate neighbors is described by the discrete Schrodinger equation 



where h is the distance between the oscillators and hZ is the lattice of points {hn : 
n — integer} . More generally, 



and the equation describing the system is given by 

iu' n (t) + A discrete u + F n (t)=0, n e hZ d , 
As h — ► 0, one obtains the continuous model. 

For F n (t) = ±\u\ 2cr u, the above equation becomes the discrete nonlinear Schrodinger 
equation. The latter is one of the prototypical differential-difference models that is 
both physically relevant and mathematically tractable. Perhaps, the most direct im- 
plementation of this equation can be identified in one-dimensional arrays of coupled 
optical waveguides [TJ |2] . These may be multi-core structures created in a slab of a 
semiconductor material (such as AlGaAs), or virtual ones, induced by a set of laser 
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1. Introduction 



iu' n (t) + h~ 2 

K(o)} e I 2 



(u n+h (t) + u n _ h (t) - 2u n (t)) + F n {t) = kehZ 




2 



A. STEFANOV AND P.G. KEVREKIDIS 



beams illuminating a photorefractive crystal. In this experimental implementation, 
there are about forty lattice sites (guiding cores), and the localized, solitary wave 
structures that this model is well known to support E] may propagate over tens 
of diffraction lengths. 

Photonic lattices [2J E] have recently provided another application of the discrete 
nonlinear Schrodinger class of models. In this case, the refractive index of a nonlinear 
medium changes periodically due to a grid of strong beams, while a weaker probe 
beam is used to monitor the localized waves. This has created a large volume of 
recent activity in the direction of understanding discrete solitons in such photonic 
lattices. 

Finally, besides its applications in nonlinear optics, the discrete nonlinear Schrodinger 
equation is a relevant model for Bose-Einstein condensates trapped in strong optical 
lattices (formed by the interference patterns of laser beams) |Hj- In this context, 
the model can be derived systematically by using the Wannier function expansions 



These applications render the discrete nonlinear Schrodinger equation a particu- 
larly relevant dynamical lattice for a variety of physical applications. 

Another example of a differential-difference equation that arises in many physical 
contexts consists of the nonlinear Klein-Gordon models. Perhaps the simplest possible 
implementation of such a lattice arises for an array of coupled torsion pendula under 
the effect of gravity [10J. However, such models are also relevant in condensed matter 
physics (e.g., describing the ffuxon dynamics in arrays of superconducting Josephson- 
junctions) [TT], as well as biophysics (e.g., describing the local denaturation of the 
DNA double strand jT5j ). It is also interesting to note that nonlinear Klein-Gordon 
models are intimately related to their Schrodinger siblings since the latter are the 
natural envelope wave reduction of the former [T2J. 

The Klein-Gordon model is of the form 



where we take h = 1 and the nonlinearity will in general be assumed to be of the 
form F n = ±\u\ 2cr u. 

In this work, we examine such Schrodinger and Klein-Gordon models as follows: 
in section 2, we give a priori energy and decay estimates for the free propagation 
in these lattices, which we subsequently prove by means of Strichartz estimates in 
section 3. In sections 4 and 5, we use these estimates to examine solutions in the 
presence of nonlinearity for the Schrodinger and Klein-Gordon cases respectively. 
As an application, we show the existence of energy thresholds for the appearance of 
localized solutions (for appropriate regimes of the parameters) and illustrate the decay 
of various lattice norms, complementing the analysis with numerical simulations. A 
number of technical details are worked out in the appendix. 



of p. 




d 2 u n (t) - A d u + u n + F n (t) = 
u n (0) = f n el 2 (Z d ), 
d t u n (0) =g n e l 2 (Z d ) 
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2. Free Evolution: Decay and Energy Estimates 
2.1. Schrodinger equation. 

Theorem 1. For the free discrete Schrodinger equation 

iu' n (t) + A d u = 
u n {0) G l\Z d ) 

one has 

(2) II {«»(*)} llf» = ||"M0)}||,a energy identity 

(3) IIK(*)}IU <c T Tj^"l|{ M n(0)}|lzi decay estimate 

Moreover, for the inhomogeneous equation 

iu' n {t) + A discrete u n + F n (t) = 
MO) e I 2 

one has the Strichartz estimates with (q,r) > 2, 1/g + dj (3r) < d/6 and 
(q, r, d) (2, oo, 3). That is 




Remark Note that the decay rate (and consequently the Strichartz estimates) is 
smaller than the usual t~ d l 2 that one has for the continuous analogue. 

In the Klein-Gordon case, similarly to the Schrodinger models, we discuss various 
relevant estimates for the free evolution. We have the following analogue of Theorem 

m 

Theorem 2. For the solutions of the one- dimensional 1 homogeneous discrete Klein- 
Gordon equation one has the a priori decay and energy estimates 

(4) ||{«t»(*)}||ja < C(||/n||/2 + \\9n\\p), energy estimate 

(5) Il{wn0)}|| /O o < ^ 1/3 ( II ,/nll/i + IbnlLO decay estimate 

1 At this stage, we have encountered some technical difficulties, when trying to extend the result 
to dimensions higher than one. In fact, the proof of the energy estimate (@J goes through in all 
dimensions, whereas the decay estimate (JSJ) boils down to the exact rate of decay of an explicit d 
dimensional integral, see While numerically, one can see that the relevant integral has the 

right rate of decay, its rigorous justification remains open. 
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Figure 1. Log-log plot of the temporal evolution of the L°° norm 
in a linear one dimensional Schrodinger lattice. The intial condition 
contains one site excited with uq = 10. The dashed line shows the best 
fit (for times 50 < t < 90) ~ r - 335 . 

For the solutions of the inhomogeneous equation, one has the Strichartz estimates 
with (q,r) > 2, l/q + l/(3r) < 1/6. That is 

111 



q/r 



ynez 



c (/(sH 



q'/f' 



dt 



3. Strichartz estimates for the Schrodinger and Klein-Gordon 

equations 

In this section, we present the proof of Theorem ^ 

Proof. (Theorem^) By rescaling in time, it suffices to consider the case h = 1. Next, 
associate {u n } < — > f(k) = Yl u n e 2mn ' k . Then, 

n£Z d 

d 

^■discrete = ^ + &j ~~ 2Jd) 

3=1 

where Sj is the shift operator in the direction of Cj. On the spaces of trigonometric 
polynomials, 

d 



A d f = J2(e 2wikj + e~ 2 ™ k > - 2)f(x) 



-4y^sin 2 (7rfej) 



/(*)■ 
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Figure 2. Log-log plot of the temporal evolution of the L°° norm in 
a linear one dimensional Klein-Gordon lattice. The intial condition 
contains one site excited with uq = 10. The dashed line shows the best 
fit (for times 20 < t < 90) ~ r - 3335 . 



Since ||K}|| /2 = II/IIl 2 ([o,i]«), & is eas y to see that U(t)u = e~ m ^^ ^f(k) is 
an isometry in I 2 . 

For the decay estimate ©, write U(t)u = e ' m ^=^^ k i) u n (0)e 2wik - n ) , whence 



(6) 



u, 



(t) = ««(0) / e-^U 



sin 2 (irkj) ^2-Ki(m—n)-k ^j^. 



Thus, (J3J), reduces to showing 



sup 

rn,n£Z a 



Ei=i sin 2 (wkj) e 2TTi(m-n)-k ^ 



0,1 



< Cmin(l, |t|- d/3 ). 



Since the integrals split into 1 D integrals, it suffices to show (after elementary change 
of variables) 

i 



sup 
seR 1 



-it(sin 2 (x)-sx) ^ 



< Cmin(l, |t 



-l/3\ 



The estimate by 1 is trivial by taking absolute values, while the estimate by |i|~ 1//3 
follows from the Van der Corput lemma, see [14] . p. 332. Indeed, we have to verify 
that the phase function tp(x) = sin 2 (a;) — sx satisfies 

max(\<p'(x)\, \ip"(x)\, \tp"'(x)\) > 1, which is verified since (ip") 2 (x) + ((p"') 2 (x) = 
4(cos 2 (2x) + sin 2 (2x)) = 4. 

Note that if s = 1, we may have max(|</?'(x)|, |<//'(x))| << 1 for x ~ 7r/4, indicating 
that this is the sharp rate of decay, at least as far as the Van der Corput lemma is 



(i 



A. STEFANOV AND P.G. KEVREKIDIS 



concerned. 

Corroborating this estimate in a direct numerical simulation, we have considered a 
linear one-dimensional lattice with h = 1. The result is shown in Fig. ^for the tem- 
poral evolution of the L°° norm. The best fit of the numerical result in the log-log plot 
of the figure very closely follows the theoretical prediction since the corresponding 
decay exponent is found to be ~ 0.335. 

The Strichartz estimates follow from the abstract Strichartz estimate by Keel and 
Tao, jT3j, which states that energy and decay estimates imply Strichartz estimates. 
The last statement is also implicit in the earlier work of Ginibre and Velo, [TSj. □ 

The proof of the Strichartz estimate for the discrete Klein-Gordon equation (DKG) 
follows closely the proof for the Schrodinger case. There are however some distinctive 
differences, which we try to highlight in the argument. 

Proof. (Theorem |2I) Let us consider the homogeneous equation first. 

The energy estimate is derived in all dimensions. Let {«„} be a solution and multiply 

both sides of the (DKG) by dtu n (t) and sum in n. We have 

d t £(<%«n) 2 - ^2(A d u n )d t u n + d t y^ j ul = 

n n n 

By the summation by parts formula, (e.g. (2.4) in ^7j, it is not hard to see that 

d 

-^2(A d u n )d t u n = <%££ \ Un ~ M n+e r | 2 / 2 - 

n r=l n 

Thus, 

d 

d t (J2( d t U n) 2 + £ £ l U " - U "+ e "! 2 + H U n) = > 
n r=l n n 

implying that 

d 

X>nW < E(^n) 2 + ^^|Wn-Wn+eJ 2 + ^^) = 

n n r=l n n 

= (E^ + EE \fn - fn + ef + £ f n ) < C(£ £ + £ £), 

n r=X n n n n 

which is (jlj). 

We have already introduced and studied the action of A disc (where for simplicity 
A disc = Ai is the discrete Laplacian in 1 D), see (JBJ). Following the same idea, 
consider {u n } < — > f = '^2 n u n e 27Tn ' k . Then, from Duhamel's formula for the wave 
equation (and also by straightforward verification), the solution to the inhomogeneous 
equation is given by 

t 

u\ (+ /i a \t a sin( V 1 - &disc) f sin((t - s)y/l - A diac ) 

u n [t) = cos^V 1 - A disc )f n H — g n + / 7= F n (s)ds. 

V 1 — A disc J V 1 — A disc 

o 
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where sin(a;), cos(a;) are expressed via the Euler's formula in terms of e %tx and a 
(continuous) function of — Adi SC is the operator given by 



(h(y/l - A disc )u) n {t) = J2 u m(°) J h(l + 4sm 2 (nk))e 



2iri(m—n)k ^j^. 



m£Z 

We have the following technical lemma, whose proof is postponed for the appendix. 

Lemma 1. The operator h{\/l — A disc ) :l p -^l p for all 1 < p < oo for all sufficiently 
smooth functions h. 

By the abtract result of Keel and Tao, the energy estimate (jH) and Lemma ^ 
matters reduce to establishing the decay estimate 

(7) 



< Ct I \\u n ||^i . 



Indeed, for Duhamel's term, we have by the above mentioned result of Keel and Tao 

t 



sm((t-sWl-A^l Fnia)da 



<||(1-A^-^ 2 F|| L , V <||F 



Vi - A disc 

LH 

where the last inequality follows from Lemma H 
As in the proof of ©, one needs to show 

i 

(8) SUp f jty/l+48m*(*k) e -2*imk dk < Cf -1/^ 

m£Z J 


which by an elementary change of variables would follow from 

l 



(9) sup / e *(0+^W- 83! )^ < Cr 1/3 , 







Thus, the phase function is tp(x) = yTT sin 2 (x) — sx and let us denote the function 
f(x) = a/1 + sin 2 (x). Denote a = arccos(2 — \/2)- 

We have that the solutions to ij)"(x) = in [0, 2n] are x = a, n — a, it + a, 2n — a. As 
we know by the Van der Corput lemma, the worst rate of decay occurs, when some 
consecutive derivatives (starting from ip') vanish at a point. In our case, we have 
that only for s = <p'(a), ip'(n — a), (p'(7r + a), (p'(2ir — a), one has ip'(x) = ip"(x) = ( 
at the points x = a, n — a, n + a, 2n — a). It is easy to check that for these values of 
s, one has max(ip' (x) , ip" (x) , ip'" (x)) > 1/2, which guarantees by the Van der Corput 
lemma a decay of at least t~ 1 ^ . 

□ 

Remark The proof of the higher dimensional analogue of such result would involve 
showing that 



(io) sup / e ^^U^)-^)) dx <cr d i\ 



seR d 

[0,1 



,1 
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Clearly, this integral does not reduce to the one dimensional case as in the Schrodinger 
case, which makes it harder object to study. Numerical simulations seem to confirm 
the validity of since in a two dimensional numerical experiment for a DKG lat- 
tice, the best fit to the decay was found to be ~ t -0 - 675 . 



(11) 



4. Applications to the discrete Schrodinger equation 
Consider 

id t u n + A d u ± \u n \ 2a u n = 



Ufc(O) G I 2 



4.1. Global solutions for small data in the regime a > 3/d. For (fTTj) . we 

establish the global existence of small solutions, provided || }||p << 1 and 
a > 3/d. 

Theorem 3. Let o > 3/d. There exists an e = e(d) > and a constant C = C(d), so 
that whenever \\{u n (0)}\\ l2 < e, a unique global solution to fjTTj) exists, which satisfies 

\\{u n (t)}\\ mr <Ce, 

for all Strichartz admissible pairs {q,r). In particular, \\u n (t)\\ lr — > as t — > oo for 
every r > 2. 

Proof. We perform an iteration procedure for the equation 

t 

u n (t) = Au n = e itA «u n (0) ± i [ e iit - s)Ad \u n \ 2r7 u n (s)ds 



in the (metric) space X = {u : sup || u ||i«;r < 2C , ||u„(0)|| i2 }, where C is the 

(q,r)— admissible 

constant in the Strichartz inequality, that is, we are seeking a fixed point of the map 
A. Clearly, by the Strichartz estimates with (q,r) and q' — l,r' — 2, we have 

sup \\Au n (t)\\ mr < C\\u n (0)\\ p +C 1 \\\u n \ 2a+1 \\ LH2 < 

(q,r)— admissible 

< C\\u n (0)\\ l2 + dw <C\\u n (0)\\ p +^11 

\\x < 

< C7|| Mn (0)|| i2 + C7 1 (2C£) 2CT+1 . 

In the above sequence of inequalities, we have used that since o > 3/d, 

(2a + l,2(2cr + 1)) is a Strichartz admissible pair and therefore controllable by the 

norm in X. Thus, for an appropriate absolute e > 0, one has 

sup \\Au n (t)\\ LlHr < 2C||« n (0)||j3, 

(q,r)— admissible 

that is A : X -> X. 
Similarly one verifies that 

||A(u n (t) — V n (t)) 1 1 _£,«/*■ rS \\ U n ~ ^7l||i2 ( T+l;2(2<T + l) ( || tt n || i 2 CT +l;2(2 ( r+l) + 1 1 V n \ \ i 2 CT + l;2(2 ( r+l) ) ^ 



^ \\u n — V n \\ L 2cr+lp(2a+l) \\ll n (0) \\p ■ 
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which by the smallness of ||w n (0) ||p implies that the map A : X —>■ X is a contraction. 
This shows the existence of a global solution of u — Au and by construction ||«|| L9/ r < 
2C\\u \\ l2 . □ 

4.2. Decay of small solutions and the Weinstein conjecture. M. Weinstein 
has proved that for a > 2/d, one has an energy excitation treshold, [T7| . i.e., that 
there exists e = s(d), so that every standing wave solution {e lAt (p n } must satisfy 
II 011 P > £■ In the same paper, he has also conjectured that for sufficiently small 
solutions, one has lim^oo ||w(t)|| ZP = for all p < oo. 

In the next theorem, we give conditions under which small solutions will actually 
decay like the free solution in the corresponding l p norms. This of course is a statement 
implying the Weinstein conjecture and is similar to a result that he has established 
for the continuous equation in an earlier paper 2 [TH] , see also the excellent expository 
paper [T5j . 

Theorem 4. Let a > 2/d and d < 2. There exists an e, so that whenever 
||-u n (0)|| Z ( 8 +2d)/(d+7) < £, one has for all p : 2 < p < (8 + 2d)/ (d + 1), 

(12) IIK(t)}lb P <C7r^- 2 )/^|| M „,(o)|| p , 

which is the generic rate of decay for the free solutions (see (fH?j) below). In particular, 
no standing wave solutions are possible under the smallness assumptions outlined 
above. 

Proof. By interpolation between (j2J) and ©, we obtain 

(13) \\{e ltAd u n (0)}\\ lP <C<t >^- 2 )/(^) \\{u n mV- 
valid for all 2 < p < oo. 

Note that from 0, we deduce that for every p > 2, ||{e i * Ad M n (0)}|| iP < £- d (p- 2 )/(3p). 

We would like to establish an a priori bound (which simultaneously implies exis- 
tence) for the solution {«„}, that is we want to place it in 

x = {{u n } : ||K(f)}IU < 2Cr^- 2 )/^)||K(0)}||^}. 

where C is the constant in the Strichartz inequality. We have by the decay estimate 

m 

\\{Au n (t)}\\ lP < C<t>- d ^- 2 ^\\u n (0)\\ lP > + 

t 

+ Ci <t _ s >( i(p-2)/(3p) \\ U n( s )\\ 2 i(2i+i)p'ds. 



2 Weinstein's result makes significant use of the pseudoconformal invariance of the continuous 
equation, which is unfortunately not present for the discrete problem. In particular, for the decay- 
ing result, he imposes the condition ||a;Mo|| i 2 < oo on the initial data. In our results, it suffices 
to assume somewhat less - in the one dimensional case: II |n| 3 / 10 M n ||, 2 << 1 and in the two 

dimensional case, one has to impose |||^| 1/ ' 2 'U n ||, 2 ± , 72 \ << b see the precise statement of Theorem 

ii in > ) 

0]for conditions in terms of the l p spaces. 
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Figure 3. Log-log plot of the temporal evolution of the / 4 (bottom 
panel), I 5 (middle panel) and I 6 (top panel) norms in a nonlinear one 
dimensional Schrodinger lattice. The intial condition contains one site 
excited with uq = 1. The dashed line shows the best fits (for times 
20 < t < 90) that are given respectively by ||w n ||j4 ~ t~ - 221 , ||tt n ||z5 ~ 

r - 257 and |K||; 6 ~r - 277 . 

Since (2a + l)p' > (4/d + 1)(8 + 2d)/(d + 7) > (8 + 2d)/(d + 1) = p for d < 2, we 
have by the inclusion V > /( 2<T+1 ) p 

\\{Au n (t)}\\ lv < c<t>- d ^- 2 ^ \\u n (0)\\ lv , + 

t 

+ Cl / <t ^ s >d( P -2)/ { 3 P ) \\ U ^ S ^ +ldS ^ 


< c < t >-^- 2 )/^) |K(o)||^ + 

+ C 2 ||m„(0)|| ;p , + J - - - - >d(p _ 2)/(3p) < - >d( p_ 2)(2ff+ i )/(3 p ) ^- 
o 

Note that d(p - 2)(2ct + l)/(3p) > d(p - 2)(4/d + l)/(3p) > 1, and therefore the 
integral term above is controlled by a constant times t~ d (p~ 2 )/( 3 f). As a consequence, 

\\{Ku n {t)}\\ lP <2C<t >-^- 2 )/^) |K(0)||^, 

provided C 3 ||u„(0)||^ < C 3 e 2r7 < C. 

The contractivity of the map A : X — > X follows in the same manner, given that 
e << 1. A fixed point argument shows that a solution exists together with the 
estimate \\{u n (t)}\\ lP < r^-2)/(3 P ). □ 

The numerical simulations performed in a d = 1, a — 3 lattice, with an initial 
condition of a single excited site with uq — 1 showed (see Fig. |3J) that the actual 
decay rate is larger than that predicted by Theorem 0] In particular, for p = 4 and 
p = 5, the theorem predicts decay rates of t~ 1//6 and t -1 / 5 respectively, while the 
numerical simulations show decay rates of t -0 - 221 and £~ - 257 respectively. For the 
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Figure 4. Same as the previous figure but for the nonlinear Klein- 
Gordon lattice. The dashed line shows the best fits (for times 20 < 
t < 90) that are given respectively by ||u„||j4 ~ r - 226 , \\u n \\{5 ~ r°- 267 
and | \u n \ |/6 ~ t~°- 292 . 

case of p — 6, the decay rate is faster and given by £ -0 - 277 (cf. with the theoretical 
prediction of t~ 2 / 9 ). 

5. Applications to the discrete Klein-Gordon equation 

For the discrete Klein-Gordon equation, as mentioned above we consider nonlin- 
earities of the form |m| 2<j m. We can deduce essentially the same results as in the case 
for the discrete nonlinear Schrodinger equation. We state them without proofs, since 
they follow by exactly the same proofs as in the Schrodinger case (by the exact same 
decay and Strichartrz estimates). 

Theorem 5. Let a > 3. There exists an e > and a constant C, so that when- 
ever ||{u n (0)}||j 2 , ||{Aw n (0)}|| Z 2 < £ there exists a unique global solution to the one- 
dimensional discrete Klein-Gordon equation with a nonlinearity \u\ 2a u, satisfying 

\\{u n (t)}\\ Lqlr <Ce, 

for all Strichartz admissible pairs (q,r). In particular, \\u n (t)\\ lr — > as t — > oo for 
every r > 2. 

In the next theorem, we establish the Weinstein excitation threshold conjecture in 
the Klein-Gordon case. 

Theorem 6. Let a > 2. There exists an e, so that whenever ||w n (0)||^/ 4 < e, 
||cVu n (0) H/5/4 < e, one has for all p : 2 < p < 5 , 

(14) IIKW}||/ P <cr^ 2 )/^)|| Mft (o)|| p ,. 

In particular, there are no small standing wave solutions to the discrete Klein-Gordon 
equation. 
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The numerical results, in this well, show a faster decay than theoreti- 

cally predicted. Furthermore, the decay is slightly faster in this model for the same 
norms in comparison with the corresponding cases for the Schrodinger equation. 
More specifically, in the Klein-Gordon model, ||wn|h 4 ~ t~ 226 , ||w n ||z 5 ~ £~ - 267 and 
||w n ||z6 ~ t -a292 . It is worth noting, however, that when we consider multiple sites 
excited by the initial condition, then we obtain decay rates which are much closer to 
the ones theoretically predicted above. More specifically, if we excite 2 nodes (rather 
than a single one), with 3 sites betwen them, then the decay rates are 
and t -0 - 270 which are considerably closer to the theoretically predicted exponents of 
1/6, 1/5 and 2/9 respectively. 

6. Appendix 

Proof. (Lemma P) we present the proof in the d dimensional case. 

We use the representation of h(y/l — Ad) to conclude that it suffices to show that for 

every N, there exists Cat, so that 



(15) 



/d 
h(l + 4^sin 2 (7r/y)e 2mm - fe dA;| < C N < m > 
m nd • J=1 



-A' 



[0,1]" 

Indeed, if we had (JTSJ), then by the expression (h(y/l — Ad)u) n = J2 m b n - m u m , we 
conclude 



h(y/l-A d ) 



u 

For the proof of (JXSJ), use integration by parts and the fact that boundary terms 

d 

disappear (due to the fact that h(l + 4 sm 2 (iikj)) is one-periodic in every variable 

i=i 

kj). We get 

„ d 

b m = — h'CL + 4 V sin 2 (7rA; 7 )) sm^irkAe^^dk. 
m jJ[o,i]« jr[ 

Clearly, one keeps integrating by parts to obtain arbitrary large rate of decay in all 
variables rrij. The lemma is proven. □ 

References 

[1] H.S. Eisenberg, R. Morandotti, Y. Silberberg, J.M. Arnold, G. Pcnnclli, and J.S. Aitchison, 
Optical discrete solitons in waveguide arrays. I. Soliton formation, J. Opt. Soc. Am. B 19 (2002) 
2938-2944. 

[2] U. Peschel, R. Morandotti, J.M. Arnold, J.S. Aitchison, H.S. Eisenberg, Y. Silberberg, T. 
Pertsch, and F. Lederer, Optical discrete solitons in waveguide arrays. 2. Dynamic properties, 
J. Opt. Soc. Am. B 19 (2002) 2637-2644. 

[3] P.G. Kevrekidis, K.O. Rasmussen, and A.R. Bishop, The discrete nonlinear Schrodinger equa- 
tion: A survey of recent results, Int. J. Mod. Phys. B 15 (2001) 2833-2900. 

[4] J.Ch. Eilbeck and M. Johansson, in Localization and Energy Transfer in Nonlinear Systems, 
L. Vazquez, R.S. MacKay, and M.P. Zorzano (eds.), (World Scientific, Singapore, 2003), p. 44. 

[5] N.K. Efremidis, S. Sears, D.N. Christodoulides, J.W. Fleischer, and M. Segev, Discrete solitons 
in photorefractive optically induced photonic lattices, Phys. Rev. E 66 (2002) 046602. 



DISCRETE SCHRODINGER/KLEIN-GORDON EQUATION 



13 



[6] A. A. Sukhorukov, Yu.S. Kivshar, H.S. Eisenberg, and Y. Silberberg, Spatial optical solitons in 

waveguide arrays, IEEE J. Quantum Elect. 39 (2003) 31-50. 
[7] F.S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smcrzi, and 

M. Inguscio, Josephson junction arrays with Bose-Einstein condensates, Science 293 (2001) 

843-846. 

[8] F.S. Cataliotti, L. Fallani, F. Ferlaino, C. Fort, P. Maddaloni, and M. Inguscio, Superfluid 
current disruption in a chain of weakly coupled Bose-Einstein condensates, New. J. Phys. 5 
(2003) 71. 

[9] G.L. Alfimov, P.G. Kevrekidis, V.V. Konotop, and M. Salerno, Wannier functions analysis of 
the nonlinear Schrodinger equation with a periodic potential, Phys. Rev. E 66 (2002) 046608. 
[10] E.N. Pelinovsky and S.K. Shavratsky, Breaking of solitary and periodic non-linear waves, Phys- 
ica 3D, 410 (1981). 

[11] A.V. Ustinov, M. Cirillo and B.A. Malomed, Fluxon dynamics in one-dimensional Josephson- 

junction arrays, Phys. Rev. B 47, 8357 (1993). 
[12] M. Peyrard and A. R. Bishop, Statistical mechanics of a nonlinear model for DNA denaturation, 

Phys. Rev. Lett., 62, 2755 (1989). 
[13] M. Peyrard and Yu.S. Kivshar, Modulational Instabilities in discrete lattices, Phys. Rev. A 46 

(1992) 3198-3205. 

[14] E. M. Stein, Harmonic analysis: Real variable methods, orthogonality, and oscillatory integrals, 

Princeton University Press, Princeton NJ, 1993. 
[15] M. Keel and T. Tao, Endpoint Strichartz Estimates, Amer. J. Math. 120 (1998), 955-980. 
[16] J. Ginibre and G. Velo, Smoothing Properties and Retarded estimates for Some Dispersive 

Evolution Equations, Comm. Math. Phys. 123 (1989), 535-573. 
[17] M.I. Weinstein, Nonlinearity 12, 673 (1999). 

[18] M.I. Weinstein, Nonlinear Schrdinger equations and sharp interpolation estimates., Comm. 

Math. Phys., 87 (1982/83), no. 4, 567-576. 
[19] M.I. Weinstein, Contemporary Math. 89, (1989). 

Atanas Stefanov, Department of Mathematics, University of Kansas, 1460 Jay- 
hawk Blvd, Lawrence, KS 66045-7523 

Panayotis Kevrekidis, Lederle Graduate Research Tower, Department of Math- 
ematics and Statistics, University of Massachusetts, Amherst, MA 01003 



